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ABSTRACT 

We present initial results from our ongoing program to image the Sunyaev-Zel'dovich (SZ) effect 
in galaxy clusters at 143 GHz using Bolocam; five clusters and one blank field are described in this 
manuscript. The images have a resolution of 58 arcsec and a radius of ~ 6 — 7 arcmin, which is 
approximately r^o — 2rsoo for these clusters. We effectively high-pass filter our data in order to 
subtract noise sourced by atmospheric fluctuations, but we are able to to obtain unbiased images of 
the clusters by deconvolving the effects of this filter. The beam-smoothed RMS is ~ 10 i^Kcmb in 
these images; with this sensitivity we are able to detect SZ signal to beyond rsoo in binned radial 
profiles. We have fit our images to beta and Nagai models, fixing spherical symmetry or allowing 
for ellipticity in the plane of the sky, and we find that the best-fit parameter values are in general 
consistent with those obtained from other X-ray and SZ data. Our data show no clear preference for 
the Nagai model or the beta model due to the limited spatial dynamic range of our images. However, 
our data show a definitive preference for elliptical models over spherical models, quantified by an F- 
ratio of ~ 20 for the two models. The weighted mean ellipticity of the five clusters is e = 0.27 ± 0.03, 
consistent with results from X-ray data. Additionally, we obtain model-independent estimates of Y500 , 
the integrated SZ y-parameter over the cluster face to a radius of r§oo, with systematics-dominated 
uncertainties of ~ 10%. Our Y500 values, which are free from the biases associated with model-derived 
Y500 values, scale with cluster mass in a way that is consistent with both self-similar predictions and 
expectations of a ~ 10% intrinsic scatter. 

Subject headings: cosmology: observation — galaxies: clusters: individual (Abell 697, Abell 1835, MS 
0015.9+1609, MS 0451.6-0305, MS 1054.4-0321) — methods: data analysis 
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1. INTRODUCTION 

Galaxy clusters are the largest collapsed objects in 
the universe, making them excellent tools for study- 
ing cosmology and the astrophysics of gravitational col- 
lapse. They are rare excursions in the matter density 
field and the formation history of clusters is closely tied 
to the composition and evolution of the universe. As 
a consequence, clusters have been used extensively to 
constrain cosmology. For example, they provided the 
first evidence that the matt er density, Q m , was insuf- 
ficient to close the universe (Bah call et all 119971 ). Ad- 
ditionally, clusters provided the most reliable informa- 
tion about the amplitude of scalar perturbations, or, 
prior to the release of WMAP resu l ts rtViana and~L iddlc 
1999t [Pieroaoli et all l200l [20031 iBorgani et al.l 120011: 
Allen et al.ll2003at iSpergel et al.ll2003D . 



Moreover, the cluster-derived measurements of cr$ can 
be combined with measurements of the normalization of 
the CMB power spectrum to constrain the tota l neutrino 
mass via the growth of density perturbations (jPierpaohl 
[200llAllen etalll2003bl: IVikhlinin et al.ll2009D . In addi- 
tion, galaxy clusters have been forming during the same 
epoch that dark energy has evolved to dominate the en- 
ergy density the universe. As a result, the number of 
galaxy clusters as a function of redshift is sensitive to 
the dark energy density, SI \ , and its equation of state, 
w (jHaiman et alj [200lt [Holder et al.l 120001 [200l . Re- 
cently, large surveys of galaxy clusters have been com- 
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pleted, or have released initial results, that constrain the 
total neutrino mas s and/or the properti e s of dark energy 
(lAllen et all 12001 IMantz et all I2010at IVikhlinin etall 
I2009t IVanderlinde et all I2010D . Additionally, clusters 
are being used to test gravity on large scales via 
studie s of their internal s t ructure and distribution in 
space (iSchmidt et~a l. 2009: Dia ferio and Ostorerol 120091: 
Martino et alj 120091: TWu et all 120101 : IMoffat and Tothl 
2010HRapetti et al.ll2010T h~ 

Galaxy clusters also provide excellent laboratories for 
studying the astrophysics of structure formation. In 
general, galaxy clusters are well-behaved objects, and 
their properties can be predicted to fairly good preci- 
sion using simple gravitat ional collapse models and self- 
similar scaling relations (Kaiser] [19860 . However, non- 
gravitational effects, such as radiative cooling, star for- 
mation, turbulence, magnetic field support, and cosmic 
ray pressure produce deviations from the simple gravita- 
tional models, and the data clearly favor models that in- 

gravita t 
INagai et al.ll2007D . 

A wide range of observational techniques are used to 
study galaxy clusters. Optical/infrared and radio mea- 
surements can be used t o study the i n dividual galaxies 
within the cluster (e.g.. iPipino et alj (|2010D : iLin et all 
( 2009): lGralla et aTT()2010l )): parameters such as velocity 
dispersion and richness can then be us ed to understand 
the global properties of the c l uster ( e . o.JMenanteau et all 
(I2010al): [Rines and DiaferTol (f20Toh; iSerra et alj (12010ft: 
iSzabo et al.l (|2010f ): lHao et al.l (|2010hT Additionally, 
X-ray observations are sensitive to the bremsstrahlung 



elude non-gravita tional processes (jKravtsov et al.1 [20061 : 
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emission from the hot gas in the intra-cluster medium 
(ICM), which co ntains ~ 90% o f the b aryonic mass of 
the cl us ter (e.g.. IVikhlinin et all (I2006D ; fG onzalc zet all 
d200l: IVikhlinin et al.l (120091): lArnaud et aLT T2009'l: 
iZhang et al.1 (|2010h ). Galaxy clusters are also efficient 
gravitational lenses, and detailed observations of the 
background galaxies can be used to determine the mat- 
ter distribution within th e cluster, ~ 90% o f which is in 
the f orm of dark matter (|Clowe et al.l 120061: lAllen et ah! 
HH). 

The ICM can also b e studied using the Sunyaey - 
Zel'dovich (SZ) effect (jSunvaev and Zel'dovichl fl972T l. 
where the background CMB photons inverse Comp- 
ton scatter off of the electrons in the ICM. Studies of 
galaxy clusters using the SZ effect are quickly matur- 
ing. The South Pole Telescope (SPT) and the Ata- 
cama Cosmology Telescope (ACT) are conducting large 
untargeted surveys and have already published cata- 
logues with dozens of c lusters (jVanderlinde et all 120101 : 
iMenanteau et a l. 2010b]); they expect to detect hundreds 
of clusters when the surveys are complete. Due to the 
redshift independence of the SZ surface brightness, the 
clusters discovered in these surveys are, on average, at 
significantly higher redshifts than those discovered with 
X-ray or optical surveys. Consequently, SZ-selected clus- 
ter catalogues are expected to play a leading role in fur- 
ther constr aining cosmological par ameters such as Q\ 
and w (e.g., iCarlstrom et all (|2002| )). 

These large-scale SZ surveys are operating based on 
the prediction that the total integrated SZ signal over 
the cluster, Y, which is proportional to the total ther- 
mal energy of the I CM, scales in a r obust way with total 
cluster mass (e.g., iKravtsov et al.l ((2006)). In practice, 
the integral does not generally extend to the edge of the 
cluster, which is effectively at > 5rsoo, and can be be 
performed in one of two ways: over a spherical volume 
using a deprojected SZ profile or over a cylindrical vol- 
ume using a projected SZ profile. Most groups have cal- 
culated the integrated SZ signal by performing the cylin- 
drical integration over the cluster face wi thin a well de- 
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though a cou ple groups have also performed the spher- 
ical i ntegral ([Mroczkowski et al.l 120091: lAndersson et all 
120101 ). With the exception of IPlagge et al.l (|2010l here^ 
after P10), these groups have exclusively used parametric 
profiles to describe the SZ signal in order to determine 
Y. The initial scaling results from these Y measurements 
are roughly consistent with the expectation that it is a 
low-scatter (~ 10%) proxy for the total cluster mass. 

Additionally, several groups have made detailed SZ 
observations of previously known galaxy clusters, and 
these data have been used to constrain the proper- 
ties of the ICM beyond r.gnn (lHalverson et al.l 120091: 
Nord et al.l 120091: iBasu et al.1 1201 0t iMroczkowski et al I 
20091: IPlagge et all l2010f )~For example, the SPT has 



measured radial profiles in 15 clusters to a significant 
fraction of the virial radius, and they find best-fit model 
parameters that are consistent with previous studies us- 
ing X-ray data P10. APEX-SZ has published detailed 
studies of three clusters, using a joint SZ/X-ray analy- 
sis to constrain mass-weighted temperature profiles be- 



yond r5op in two of these clusters ( Halverson et al.ll2009l : 
iNord et all 120091: IBasu et all 120101 ) . Note that, in addi- 
tion to providing complementary information about the 
ICM in these previously studied clusters, these SZ data 
will help constrain and improve systematic uncertainties 
in blind SZ surveys (e.g., constraining the typical SZ 
profiles and resulting detection bi ases due to the cool- 
roperties of the central gas (|Pipino and Pierpaolil 
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High resolution (~ 10 arcsec) 
ments are also now b e ing made 
TANG dMason et all 120101: iKo-rngut et al.l 120101 ) and 
SZA/CARMA (jC^xlstrorn_elalJ|2002j). These observa- 
tions will provide insights on the internal structure of 
clusters that will be complementary to the information 
obtained via X-ray observations. Among other appli- 
cations, these data can be used, in combination with 
optical observations, to infer the merging history and 
evolution of the cluster ()Croston et al.|[2008() . 

In this work, we present Bolocam SZ observations of 
five massive galaxy clusters. We are able to measure SZ 
signal in our two-dimensional images to ~ rsoo and to 
well beyond r5oo in azimuthally averaged radial profiles. 
Using these data, we are able to constrain the broad mor- 
phologies of these clusters and compute observables such 
as Ygz in a model-inde pendent way. A comp anion pa- 
per to this manuscript, lAmeglio et all ([2010), presents 
joint X-ray/SZ deprojections of density and temperature 
profiles for the se same five clusters using Bolocam and 
Chandra data (jAmeglio et alJl2007l ). 

2. CLUSTER OBSERVATIONS WITH BOLOCAM 

Bolocam is a mm-wave imaging camera that operates 
from the Caltech Submillimeter Observatory (CSO) with 
144 bolometric det ectors covering a circular 8 arcmi n 
field of view (FOV) (|Glenn et al.lll998t IHaig et~aIll20T)l . 
For the observations described in this manuscript, Bolo- 
cam was configured to observe at 143 GHz. To image 
the c lusters, we scanned the telescope in a Lissajous pat- 
tern (jKovacs et al.ll2006l ). where the telescope is driven in 
two orthogonal directions using sine waves with incom- 
mensurate periods. We used an amplitude of 4 arcmin for 
the sinusoids, which were oriented along the RA and dec 
directions with periods of 6.28 and 8.89 sec. These pa- 
rameters were chosen to keep the FOV on the cluster cen- 
ter 100% of the time while scanning as fast as possible at 
the CSO (~ 4 arcmin/sec) to modulate the cluster signal 
above the low-frequency atmospheric noise. An example 
of the resulting integration time per pixel is given in Fig- 
ure [TJ The data were collected via 10-minute-long obser- 
vations, with the periods of the two orthogonal sinusoids 
exchanged between observations. Typically, we complete 
~ 100 observations per cluster, which corresponds to ap- 
proximately 50 ksec and yields a beam-smoothed RMS 
of ~ 10 ijKcmb (see Table [1]). 

In this paper, we present the results from five clusters 
and one blank field, which are described below. Due to 
the size of our resulting maps (r ~ 6 — 7 arcmin), we 
have chosen to focus primarily on high-redshift clusters, 
which have virial radii within the extent of our maps. All 
of the clusters presented in this manuscript are beyond 
z = 0.25, and three of the five are beyond z = 0.50. For 
reference, t^qq lies within the extent of our map for a 
typical 10 15 M Q cluster at z = 0.25, and at z = 0.50, the 
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TABLE 1 

Cluster, properties 



target 


RA 


dec 


redshift 


Bolocam time (ksec) 


RMS (vKcmb) 


r 500 (Mpc) 


M gaSt500 (M ) 


Abell 697 


08:42:58 


+36:21:56 


0.28 


52 


8.9 


1.65 + 0.09 


19.6 + 2.7 x 10 13 


Abell 1835 


14:01:02 


+02:52:42 


0.25 


50 


8.7 


1.49 + 0.06 


14.1 ± 1.2 x 10 13 


MS 0015.9+1609 


00:18:34 


+16:26:13 


0.54 


38 


10.2 


1.28 + 0.08 


17.5 + 1.9 x 10 13 


MS 0451.6-0305 


04:54:11 


-03:00:53 


0.55 


53 


7.7 


1.45 + 0.12 


15.6 + 2.2 x 10 13 


MS 1054.4-0321 


10:56:59 


-03:37:34 


0.83 


66 


6.7 


1.07 + 0.13 


11.5 + 2.4 x 10 13 


SDS1 


02:18:00 


-05:00:00 




37 


9.1 







Note. — A list of the clusters presented in this manuscript. From left to right the columns give the RA and dec of the cluster 
in J2000 coordinates, the redshift of the cluster, the amount of Bolocam integration time, the median RMS per beam-smoothed 
pixel in the Bolocam map , the radius of the cluster, and the mass of the c luster . The values for rsoo and M gas 500 were taken 
from IMantz et al.l Q2010bfl (Abell 697 and Abell 1835) and lEttori eTaT] p009h (MS 0015.9+1609, MS 0451.6-0305, and MS 
1054.4-0321) 
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Fig. 1. — Integration time per pixel, relative to the maximum 
integration time, for MS 0451.6-0305. Our model fits include all of 
the data within a circular region with a minimum integration time 
of 25% of the peak integration time, which corresponds to 6 — 7 
arcmin in radius. The red box, with 10 arcmin sides, denotes the 
region used for deconvolution of the processing transfer function. 
The minimum relative integration time within this region is also 
> 25%. 



virial radius lies within our map. Note that, throughout 
this work, we assume a ACDM cosmology with Q m = 0.3, 
Q A = 0.7, and H = 70 km/s/Mpc. 

Abell 697: Abell 697 is a cluster undergo ing a complex 
merg er event along the line of sight (jGirardi et al.l 

Abell 1835: Abell 1835 is a relaxed cluste r with 
a strong cooling flow (jPeterson et al.l 120011 : 
ISchmidt et al.ll200lh . 

MS 0015.9+1609: MS 0015.9+1609 is a triaxial clus- 
ter that is elongated along the line of sight and 
has an anomalously h igh gas mass fraction of 27% 
(jPiffaretti et al.ll200l . 

MS 0451.6-0305: MS 0451.6-0305 is a cluster that 
is not quite in gravitational eq uilibrium, with 
a slig htly elongated X-ray profile (|Donahue et al.l 
2003). 



MS 1054.4-0321: MS 1054.4-0321 is a cluster under- 
going a merger, as evidenced by the presence 
of two distinct sub-clumps in the X-ray image 
(jJeltema et al.ll200lt Uee et al.|[2005l ). 

SDS1: Our SDS1 map is centered in the middle of the 
Subaru/XMM Deep Survey (SXDS) field. The 
deep XMM-Newton survey of this field reveals only 
3 clusters within our map, all near the edge, and 
the largest of w hich has a virial mass o f M200 = 
0.8 x 10 13 M (jFinoguenov et alJl20l"ch . There- 
fore, SDS1 is approximately free of signal from the 
SZ effect. 

3. DATA REDUCTION 

In general, our data red u ction followed the procedure 
described in iSayers et al.1 (|2009L hereafter S09), with 
some minor modifications. We briefly describe the tech- 
niques below, along with the changes relative to S09. 

3.1. Calibration 

Bright quasars located near the clusters were observed 
for 10 minutes once every ~ 90 minutes in order to 
determine the offset of our focal plane relative to the 
telescope pointing coordinates. These observations were 
used to construct a model of the pointing offset as a func- 
tion of local coordinates (az,el), with a single model for 
each cluster. The uncertainty in the pointing models is 
< 5 arcsec. This pointing uncertainty is quasi-negligible 
for Bolocam's 58 arcsec FWHM beams, especially for 
extended objects such as clusters. We made two 20- 
minute-long o bservations eac h night of Uranus, Neptune, 
or a source m Sandcll (1994) for flux calibration. Using 
the quiescent detector resistance as a proxy for detec- 
tor responsivity and atmospheric transmission, we then 
fit a single flux-calibration curve to the entire data set. 
We estimate the uncertainty in our flux calibration to be 
4.3%, with the following breakdown: 1.7% from the Rudy 
temper ature model of Mars scaled to measured WMAP 
value s (lHalverson et a,l.ll2009t lWrightill976t iGriffin et all 
19861: lRudvetal.1 119871: iMuhleman and Bergd 119911: 
Hill et all |2009D. 1 5% in the Uranus /Nept une model ref- 
erenced to Mars (jGriffin and Ortonl fl993] ) . 1.4% due to 
variations in atmospheric opacity (S09), 3.1% due to un- 
certainties in the solid angle of our point-spread function 
(PSF) (S09), and 1.5% due to measurement uncertain- 
ties (S09). 

3.2. Atmospheric noise subtraction 
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Fig. 2. — Timestream noise PSD for a typical Bolocam detec- 
tor. The black curve shows the raw PSD recorded by the detector; 
spectral lines at the fundamental scan frequencies are clearly seen 
above the broadband atmospheric noise. The red curve shows the 
noise PSD after subtracting the atmospheric noise using the aver- 
age signal over the FOV. This timestream is then high-pass filtered 
at 250 mHz to produce the green PSD. Note that there is very little 
cluster signal above ~ 2 Hz, where there are some spectral lines 
due to the readout electronics. The dashed horizontal line provides 
an estimate of the photon, or BLIP, noise. 

The raw Bolocam timestreams are dominated by noise 
sourced by fluctuations in the water vapor in the atmo- 
sphere, which have a power spectrum that rises sharply 
at low frequencies. In order to optimally subtract the 
atmospheric noise, we have used a slightly modified ver- 
sion of the average subtraction algorithm described in 
ISavers et"al~l (|2010[ hereafter S10). We have modified 
the S10 algorithm because these cluster data contain ad- 
ditional atmospheric noise caused by the Lissajous scan 
pattern. Since we are scanning the telescope parallel to 
RA and dec, the airmass we are looking through is con- 
stantly changing. As a result, our data contain a large 
amount of atmospheric signal in narrow bands centered 
on the two fundamental scan frequencies. 

Following the algorithm in S10, we first create a tem- 
plate of the atmosphere by averaging the signal from all 
of our detectors at each time sample (i.e., the average 
signal over the FOV). In S10, this template is subtracted 
from each detector's timestream after weighting it by the 
relative gain of that detector, which is determined from 
the correlation coefficient between the timestream and 
the template. We use a single correlation coefficient for 
each detector for each 10-minute-long observation. How- 
ever, a significant fraction of the atmospheric noise at the 
fundamental scan frequencies remains in the data after 
application of the S10 algorithm, indicating that we have 
slightly misestimated the correlation coefficients. There- 
fore, we modified the S10 algorithm to compute the cor- 
relation coefficients for the template based only on the 
data within a narrow band centered on the two funda- 
mental scan frequencies. The atmospheric noise power 
in these narrow frequency bands is roughly an order of 
magnitude above the broadband atmospheric noise at 
nearby frequencies; consequently, the data in these nar- 
row bands provide a high signal-to-noise estimate of each 
detector's response to atmospheric signal. The narrow- 
band atmospheric noise features are completely removed 
using this modified S10 algorithm, and the amount of 
residual broadband atmospheric noise is slightly reduced 
compared to the results from the original S10 algorithm. 

After applying this average subtraction algorithm to 



the timestream data, we then high-pass filter the data 
according to 

1 

T=l- 



1 + (10//A-1)k 

with /o = 250 mHz and k = 8. The value of fo was cho- 
sen to maximize the spatially-extended S /N for the typ- 
ical cluster in our sample based on tests with fo varying 
from to 400 mHz, and the value of k was chosen to pro- 
duce a sharp cutoff with minimal ringing. Figure [2] shows 
a typical pre and post-subtraction timestream noise PSD. 

3.3. Transfer function of the atmospheric noise filtering 

In addition to subtracting atmospheric noise, the FOV- 
average subtraction and timestream high-pass filter also 
remove some cluster signal. Since we use the data 
timestreams to both determine the atmospheric fluctua- 
tion template and the correlation coefficient of each de- 
tector's data timestream with the template, the FOV- 
average subtraction acts on the data in a non-linear way. 
Consequently, its impact on the data depends on the clus- 
ter signal. As described below, we quantify the effects of 
the FOV-average subtraction and the timestream high- 
pass filter via simulation by processing a known cluster 
image through our data-reduction pipeline. These sim- 
ulations are computation-time intensive, and we find, in 
practice (see Section PO]) . that the filtering is only mildly 
dependent on the cluster signal. Thus, in the end, we 
determine the effects of the filtering for a particular clus- 
ter's data set using the cluster model that best fits those 
data. We use the term transfer function to describe the 
effect of the filtering, although this terminology is not 
rigorously correct because the filtering depends in the 
cluster signal. 

To compute the transfer function, we first insert a 
simulated, beam-smoothed cluster profile into our data 
timestreams by reverse mapping it using our pointing in- 
formation. These data are then processed in an identical 
way to the original data, and an output image, or map, 
is produced. When processing the data-plus-simulated- 
cluster timestreams we use the FOV-average subtraction 
correlation coefficients that were determined for the orig- 
inal data. This ensures that the simulated cluster is pro- 
cessed in an identical way to the real cluster in our data. 
In the limit that the best-fit cluster model is an accurate 
description of the data, this process is rigorously correct. 
The original data map is then subtracted from this data- 
plus-simulated-cluster map to produce a noise-free image 
of the processed cluster. In Figured we show an example 
cluster image, along with the noise-free processed image 
of the same cluster. The Fourier transform of this pro- 
cessed cluster image is divided by the Fourier transform 
of the input cluster to determine how the cluster is fil- 
tered as a function of 2-dimensional Fourier mode (i.e., 
what we term the transfer function, see Figure 2]). 

At small angular scales, there is very little signal in the 
beam-smoothed input cluster, and numerical noise pre- 
vents us from accurately characterizing the transfer func- 
tion at these scales. The transfer function is expected to 
be unity at small angular scales, and is asymptoting to 
this value at the larger scales where we can accurately 
characterize it. Therefore, we set the transfer function 
to a value of 1 for u > 0.75 ajcmin . Deconvolving the 
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Fig. 3. — Images of the best-fit spherical Nagai model for MS 0451.6-0305. The left image is the model and the right image is the model 
after being processed through our data reduction pipeline, which high-pass filters the image in a complex way. This filtering significantly 
reduces the peak decrement of the cluster and creates a ring of positive flux at r > 2 arcmin. Note that the processed image is not quite 
azimuthally symmetric. 
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Fig. 4. — The magnitude of the transfer function for MS 0451.6- 
0305 as a function of Fourier wavenumber u = 1/A. At large 
scales, or small u, the measurement error is negligible and the 
error bars provide an indication of the azimuthal variation. At 
u > 0.75 arcmin , the measurement error becomes non-negligible 
and we set the transfer function equal to 1. Note that this az- 
imuthally averaged transfer function is for display purposes only; 
we have used the full two-dimensional transfer function throughout 
our analysis. 

processed cluster image using the transfer function (see 
Section [5]) , which has been approximated as 1 at small 
angular scales, produces an image that is slightly biased 
compared to the input cluster. The residuals between 
these two images are approximately white, with an RMS 
of < 0.1 fiKcMB- This transfer- function-induced bias is 
negligible compared to our noise, which has an RMS of 

~ 10 flKcMB- 

As noted above, the transfer function (weakly) de- 
pends on the profile of the cluster; larger clusters are 
more heavily filtered than smaller clusters. Therefore, 
we determine a unique transfer function for each cluster 
usin g the best-fit ellip tical Nagai model for that clus- 
ter ijNagai et all 120071 hereafter N07). The details of 
this fit are given in Section Since the transfer func- 
tion depends on the best-fit model, and vice versa, we 
determine the best-fit model and transfer function in an 
iterative way. Starting with a generic cluster profile, we 



first determine a transfer function, and then fit an ellip- 
tical Nagai model using this transfer function (i.e., the 
Nagai model parameters are varied while the transfer 
function is held fixed). This process is repeated, using 
the best-fit model from the previous iteration to calculate 
the transfer function, until the best-fit model parameters 
stabilize. This process converges fairly quickly, usually 
after a single iteration for the clusters in our sample. 
The model dependence of the resulting transfer function 
is quantified in Section ROI 

3.4. Noise estimation 

In order to accurately characterize the sensitivity of 
our images, we compute our map-space noise directly 
from the data via 1000 jackknife realizations of our clus- 
ter images. In each realization, random subsets of half 
of the ~ 100 observations are multiplied by —1 prior to 
adding them into the map. Each jackknife preserves the 
noise properties of the map while removing all of the as- 
tronomical signal, along with any possible fixed-pattern 
or scan-synchronous noise due to the telescope scanning 
motion^. Since these jackknife realizations remove all as- 
tronomical signal, we estimate the amount of astronom- 
ical noise in our images separately, as described below. 
After normalizing the noise estimate of each map pixel 
in each jackknife by the square root of the integration 
time in that pixel, we construct a sensitivity histogram, 
in iiKcmb-s 1 ^ 2 , from the ensemble of map pixels in all 
1000 jackknifes. The width of this histogram provides an 
accurate estimate of our map-space sensitivity (see Fig- 
ure [5]). We then assume that the noise covariance matrix 
is diagonal and divide by the square root of the integra- 
tion time in each map pixel to determine the noise RMS 
in that pixel. This method is analogous to the one used 
in S09, where it is described in more detail. 

There is a non-negligible amount of noise in our 

4 We show that there is no measurable fixed-pattern or scan- 
synchronous noise in our data later in this section and in Sec- 
tion U£2] 

5 T his ap proximation is justified for our processed data maps in 
Scction l6.2l Note that the approximation fails for our deconvolved 
images (sec Section , which contain a non- negligible amount of 
correlated noise. We describe how this correlated noise is accounted 
for in our results in Section [5] 
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Fig. 5. — Histogram of the per-pixel sensitivity for our maps of 
the blank-field SDSl. The black line shows the average histogram 
for each of the 1000 jackknife realizations of the data, with noise 
from th e CM B and unresolved point sources added as described in 
Section 13.41 The sensitivities derived from the noise realizations 
are well described by a Gaussian fit, with the fit quality quanti- 
fied by a PTE of 0.44. Overlaid in red is the histogram for our 
map of SDSl, which is also well described by the Gaussian fit 
to the noise realizations (PTE = 0.63), indicating that our noise 
model adequately describes the data. The dashed blue line shows 
the best-fit Gaussian to the noise realizations and has an RMS of 
4.02 hKcmb-s 1 ^. 

maps from two main types of astronomical sources: 
anisotropies in the CMB and unresolved point sources. 
The South Pole Telescope (SPT) has recently published 
power spectra for both of these sources at 150 GHz over 
the range of angular scales probed by our maps, so we use 
their measur ements to estimate the astronomica l noise 
in our maps (ILueker et all [Mot IHall et al.l[2010h . Note 
that these measurements inherently include all astro- 
physical effects, such as lensing and clustering of point 
sources. Using the SPT power spectra, we generate simu- 
lated maps of our cluster fields assuming Gaussian fluctu- 
ations. Note that this is a poor assumption for both the 
SZ-sourced CMB fluctuations and the signal sourced by 
background galaxies that are lensed by the cluster. How- 
ever, since the noise power from both of these sources is 
quasi-negligible compared to the total noise in our im- 
ages, any failure of our Gaussian assumption will have 
a minimal impact on our results. The good match of 
the SDSl data to our noise model validates our Gaus- 
sian assumption (see Figure [5]). These simulated images 
of the CMB plus point sources are then reverse mapped 
into our timestream data and processed to estimate how 
they will appear in our cluster maps. We then add these 
processed astronomical realizations to our jackknife real- 
izations to provide a complete estimate of the noise in our 
images. Sensitivity estimates from our signal-free map of 
SDSl agree well with the sensitivity estimates from our 
noise maps, indicating that there are no additional noise 
sources that have not been included in our estimate (see 
Figure[5]). The properties of this noise are described in 
more detail in Section 16.21 In general, the noise from 
astronomical sources is quasi-negligible compared to the 
other noise in our images. However, the large-scale corre- 
lations from CMB fluctuations are non-negligible in our 
deconvolved images (see Table E|) . 

Note that we have chosen to model the signal from 
point sources as an additional noise term in our maps 
rather than attempting to subtract any individual point 
sources. Part of the motivation for this approach is the 
fact that we do not detect any point sources in our cluster 



or blank-field images. For reference, P10 mapped 15 clus- 
ters to a similar depth using the SPT and only detected 
two point sources; the combined area of the SPT images 
is approximately 20 times the combined area of our 6 im- 
ages. Although there are several known SMGs and AGNs 
in our images (|Zemcov et al.1 120071 llvison et al.1 120001 : 
ICoorav et al.lll998D . even the brightest sources will have 
a flux of ~ 10 /iKcMB in our observing band, rendering 
them undetectable given our noise RMS of ~ 10 [iKcmb- 
Additionally, it is not possible to reliably estimate the 
flux of these sources in our observing band since most 
have only been detected at one or two wavelengths, gen- 
erally separated from our observing band by more than 
a factor of 2 in wavelength. 

4. CLUSTER MODEL FITTING 

4.1. The SZ effect 

As mentioned above, the SZ effect involves CMB pho- 
tons inverse Compton scattering off of hot electrons in 
the ICM. Since the electrons are many orders of magni- 
tude hotter than the CMB photons, there is, on average, 
a net increase in the energy of the photons. The classi- 
cal distortion relative to the blackbody spectrum of the 
CMB is given by 

*/ ^ eX + 1 a 

where x = hv /ksTcMB- The magnitude of the distor- 
tion is proportional to the product of the density and 
temperature of the electrons in the ICM projected along 
the line of sight. The frequency-dependent temperature 
change of the CMB is given by 

Tsz = f(x)yT C MB, 

where 



and n e and T e are the density and temperature of the 
ICM electrons. Relativistic correction s can be include d 
by multiplying f(x) by (1 + 5(x,T e )) (jltoh et al.1 1 1991 . 
Note that the X-ray brightness of the ICM is propor- 
tional to rigT e ; the differing density and tempera- 
ture dependence of the SZ and X-ray signals makes 
them highly compleme ntary probes of the ICM (e.g., 
IBonamente et all (|2006h . hereafter B06). 

4.2. Models 

We have fit our d ata to two types of models: an isother- 
mal beta model ( Cavalie re and Fusco-Femianol 119761 
[1971 and the pressure profile proposed by N07, here- 
after the Nagai profile. The isothermal beta model has 
been used extensively to describe X-ray and SZ mea- 
surements of the ICM, and our beta model fit parame- 
ters can be directly compared to these previous results. 
However, the beta model provides a poor description of 
deep X-ray data, which gen erally have a cuspier co re and 
stee per outer profi l e (e.g ., iVikhlinin et al.l ((2006), N07, 
and lArnaud et al.1 ()2009t )V In contrast, the Nagai pres- 
sure profile, which is a generalization of the NFW dark 
matter profile, is able to describe deep X-ray observatio ns 
over a wide range of scales (N07, lArnaud et all ([2009)). 
Therefore, we have also fit our data to the Nagai model. 
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In practice, the beta and Nagai models are highly degen- 
erate over the range of angular scales to which our data 
are sensitive (1 — 12 arcmin). 
Specifically, the isothermal beta model is described by 



P = 



Po 



where p is the pressure profile, po is the pressure normal- 
ization, r c is the core radius, and /3 is the power law slope. 
The beta model can be analytically integrated along the 
line of sight to give the observed SZ signal, with 



Tsz — 



f(x)ypT C MB 

(1 + r 2/ r 2)(3/3-l)/2 



■ST, 



where Tsz is the SZ signal in our map (in plKcmb) and 
?/o is the central Comptonization. We fix the value of (3 
at 0.86 for all of our fits; this is the best-fit value found in 
P10 for SZ datsQ. Since our data timestreams are high- 
pass filtered, we are not sensitive to the DC signal level 
in our images. Note that because the timestream data, 
rather than the map, are high-pass filtered, the DC signal 
level of the map is in general not equal to 0. However, 
the DC signal of the map is not physically meaningful 
due to the filtering and must therefore be included as 
a free parameter, ST, in the model fits. We have also 
generalized the beta model to be elliptical in the plane 
of the sky with 



Tsz — 



f(x)yoT ( 



CMB 



+ £ 



(1- 



ST, 



where n is oriented along the major axis, described by 
a position angle of 9 (in degrees east of north), r 2 is 
orthogonal to the major axis, and e is the ellipticity. 
The Nagai model is described by 



P 



Po 



(r/r s ) c [l + {r/r s f] {B - C)/A ' 



(1) 



where p is the pressure, po is the pressure normaliza- 
tion , r s is the scale r adius (r s = t^qq/c, with c ~ 
1.2 (jArnaud et alj |2009)). and A, B, and C are the power 
law slopes at intermediate, large, and small radii com- 
pared to r s . We have fixed the values of A, B, and 
C to the best-fit values found in lArnaud et al.1 (2009) 
(1.05,5.49,0.31). The Nagai model is not analytically in- 
tegrable, so we numerically integrate p along the line 
of sight to determine Tsz- We have also generalized 
the Nagai model to be elliptical in the plane of the 
sky, using the same notation as our elliptical general- 
ization of the beta model. Note that, in all of our fits, 
we have corrected for relativi s tic eff ects via the approx- 
imations given in lltoh et all JI998) using gas temper- 
ature es_dmate^_from Xnray observations of thes e clus- 
ters ((C^gnoioiFin [Hi EeEe^eFaDEOOil). The 
relativistic corrections are ~ 5% for the typical temper- 
atures in these massive clusters (~ 10 keV). 

We find that both the beta model and the Nagai model 
adequately describe our data, with the exception of Abell 

6 As P10 point out, this value for /3 is larger than those 
gener ally found from X- ray data, in agreement with simula- 
tions (Hallman et al. 2007). Additionally, it suggests that the ICM 
temperature is falling with increasing radius, in agreement with 
simulations and data (e.g., N07). 



697, and neither model is preferred with any significance. 
However, the elliptical models (x 2 / D OF = 5605/5413) 
provide a much better fit to the full ensemble of our 
cluster data compared to the spherical models (% 2 /DOF 
= 5813/5423, see Tables H and [3j . The F-ratio for these 
two fits is 20.1 for = 10 and v-i = 5413 degrees 
of freedom (|Bevington a nd Robinson 1992]), correspond- 
ing to a probability of < 10 -36 that we would obtain 
data with equal or greater preference for elliptical mod- 
els if the clusters were spherical. If we neglect Abell 
697, which is not well described by either model, then 
the F-ratio is 12.2 with a corresponding probability of 
< 10~ 16 . The weighted mean ellipticity of the five clus- 
ters is e = 0.2 7 ± 0.03, consistent wi t h results from X- 
ray da ta (e.g., Maugha n et al.l (|2008f ); I5e Filippis et al.l 
(2005)). Our inability to distinguish between the Nagai 
and beta models is likely due to the limited spatial dy- 
namic range of our images (~ 1 — 12 arcmin). For the 
clusters in our sample, we are insensitive to the core, 
r < 0.15r 50 o — 0.5 arcmin, and the outskirts of the clus- 
ter r > 2r5oo — 10 arcmin, which means we are only 
sensitive to a single power law exponent in the Nagai 
model, A. For this reason, the Nagai model is highly 
degenerate with the beta model for the angular scales 
probed by our data. 

4.3. Fitting procedure 

Our cluster images are filtered by both the atmospheric 
noise subtraction and the Bolocam point-spread func- 
tion. Therefore, prior to fitting a model to our data, 
we need to filter the cluster model in an identical way. 
First, we generate a 2-dimensional image using the clus- 
ter model, computed directly from the isothermal beta 
model and via a line-of-sight integration of the Nagai 
pressure model. Next, the model image is convolved 
with the Bolocam point-spread function and the mea- 
sured transfer function. In practice, the transfer function 
convolution is performed via multiplication in Fourier 
space. This filtered model is then compared to our data 
map, using all of the map pixels contained within a ra- 
dius where the minimum coverage is greater than 25% of 
the peak coverage. This radius is typically between 6 and 
7 arcmin. We use an iterative least-squares technique to 
determine the best-fit parameters for each model; we esti- 
mate the uncertainty on each parameter via the standard 
deviation of the best-fit parameter values estimated from 
noise realizations with model clusters added to them (see 
Section |6~2"|). 

For each model, we fit both a scale radius (r c for the 
beta model and r s for the Nagai model) and a normaliza- 
tion (uq for the beta model and po for the Nagai model) . 
Additionally, as described above, we fit for the observa- 
tionally unconstrained DC signal level of the map, ST. 
We also fit for a centroid offset relative to the X-ray 
pointing center. The offsets for the five clusters range 
from ~ — 20 ± 5 arcsec, indicating there are no major 
differences in the X-ray and SZ centroids. Finally, when 
we allow the model to be elliptical in the image plane, we 
fit for the ellipticity e and position angle 9. A complete 
list of the best-fit parameters for spherical and elliptical 
versions of both models is given in Tables [2] and [3l 

4.4. Discussion and comparison to previous results 



TABLE 2 
Nagai model fit parameters 



cluster po 


/ 10 -il_erg , 

\ cm 3 ' 


r s (arcmin) 


C500 






(deg) 


X 2 /DOF 


PTE X 2 


PTE ■ 

r L ^sim 


elliptical Nagai model 




















Abell 697 


9.3 + 1.3 


6.9 + 1.0 


0.93 ± 


14 


0.37 + 0.05 


-24 + 4 


1289/1117 


0.00 


0.00 


Abell 1835 


8.1 ± 1.7 


6.7+1.5 


0.94 ± 


21 


0.27 + 0.07 


-16 ± 10 


966/945 


0.31 


0.22 


MS 0015.9+1609 


6.7 + 1.4 


5.5 + 1.1 


0.62 + 


13 


0.24 ± 0.08 


68 ± 12 


1079/1117 


0.79 


0.73 


MS 0451.6-0305 


8.7 + 1.6 


4.7 + 0.7 


0.80 ± 


11 


0.26 + 0.06 


85 + 7 


1188/1117 


0.07 


0.08 


MS 1054.4-0321 


6.5 + 1.9 


3.6 + 1.0 


0.64 ± 


19 


0.09 + 0.07 


-1 ± 32 


1084/1117 


0.75 


0.72 


spherical Nagai model 




















Abell 697 


11.0 + 1.8 


4.6 + 0.8 


1.40 + 


22 






1399/1119 


0.00 


0.00 


Abell 1835 


8.9 + 2.1 


5.1 ± 1.2 


1.24 + 


27 






1007/947 


0.08 


0.07 


MS 0015.9+1609 


6.0 + 1.1 


5.4 + 1.1 


0.63 ± 


12 






1100/1119 


0.66 


0.61 


MS 0451.6-0305 


10.6 + 2.1 


3.5 + 0.5 


1.08 + 


15 






1220/1119 


0.02 


0.02 


MS 1054.4-0321 


6.0 + 1.8 


3.6 + 1.0 


0.64 ± 


16 






1087/1119 


0.75 


0.75 



Note. — Table of the best- fit parameters and la uncertainties for our Nagai model fits with the power-law exponents fixed 
to the best-fit values found in Arnaud ct al. (2009) (1.05,5.49,0.31). From left to right the columns give the normalization of 
the Nagai profile, po, the scale radius of the major axis, r 3 , the concentration parameter, C500, the ellipticity, e, the position 
angle of the major axis in degrees east of north, 8, the x an d DOF for the fit, and the goodness of fit quantified by the 
probability to exceed the given x 2 /DOF based on the standard x 2 probability distribution function and also empirically by the 
fraction of our 1000 noise realizations that produce a larger x 2 value when a model cluster is added to t hem (see Sectio n 16.211 . 
The values of C500 for the spherical fits can be co mpar ed to the nominal value of 1.17 that is given in Arnaud ct al. (2009) 
based on spherical fits of X-ray data. See Section 16.21 for a description of how the parameter uncertainties are calculated. 
Note that we have included the 4.3% uncertainty in our flux calibration in the error estimates for pp. 



TABLE 3 
Beta model fit parameters 



cluster 



yo (10 4 ) r c (arcsec) r c /r 500 



6 (deg) 



X 2 /DOF PTE 2 PTE S 



elliptical beta model 



Abell 697 


2 


93 + 


26 


98 ± 


10 





26 + 


03 





38 ± 0.04 


-10 + 2 


1288/1117 





00 





00 


Abell 1835 


2 


49 + 


31 


79 + 


11 





21 + 


03 





28 + 0.07 


-9 + 10 


970/945 





28 





20 


MS 0015.9+1609 


2 


59 + 


27 


85 ± 


12 





43 + 


07 





24 + 0.07 


70 + 12 


1078/1117 





79 





73 


MS 0451.6-0305 


2 


89 + 


22 


68 ± 


8 





31 + 


01 





26 ± 0.06 


80 + 8 


1193/1117 





06 





06 


MS 1054.4-0321 


2 


13 + 


22 


58 ± 


12 





41 + 


10 





09 ± 0.07 


7 + 33 


1083/1117 





76 





73 


spherical beta model 
































Abell 697 


2 


74 + 


21 


72 + 


8 





19 + 


02 








1397/1119 





00 





00 


Abell 1835 


2 


46 + 


31 


63 ± 


8 





17 + 


02 








1008/947 





08 





08 


MS 0015.9+1609 


2 


54 + 


29 


80 + 


11 





40 + 


05 








1099/1119 





66 





61 


MS 0451.6-0305 


2 


92 + 


24 


53 + 


7 





24 + 


03 








1222/1119 





02 





02 


MS 1054.4-0321 


2 


13 + 


22 


56 + 


10 





40 + 


07 








1085/1119 





76 





76 



Note. — Table of the best-fit parameters and lcr uncertainties for our beta model fits with the power law exponent 
/9 set to the best-fit value found in P10 (0.86). From left to right the columns give the normalization of the beta profile, 
yo, the core radius of the major axis, r Ct the relative value of the core radius, r c /r5oo, the ellipticity, e, the position 
angle of the major axis in degrees east of north, 9, the x 2 an d DOF for the fit, and the goodness of fit quantified by the 
probability to exceed the given x 2 /DOF based on the standard x 2 probability distribution function and also empirically 
by the fract ion of our 1000 noise realizations that produce a larger x 2 value when a model cluster is added to them (see 
Section l6.2H . The values of r c /r^oo for the sphe rica l fits can be compared to the nominal value of 0.20 given in P10 based 
on spherical fits of SPT SZ data. See Section 16.21 for a description of how the parameter uncertainties are calculated. 
Note that we have included the 4.3% uncertainty in our flux calibration in the error estimates for yo. 



All of these clusters have been studied extensively at 
a wide range of wavelengths, providing us with a large 
number of published results to compare our model fit 
parameters to. In general, our results agree well with 
those found from previous studies. In particular, when 
we fit a spherical beta model to our data using the values 
of /3 given in B06, which were determined using Chandra 
X-ray data and 30 GHz interferometric SZ data, our 
best-fit values for yo agree quite well, with the exception 
of Abell 697 (see Table EJ). The best-fit values we find 
for r c are also consistent with those determined in B06, 
with the exception of Abell 1835 and MS 0451.6-0305. 
Additionally, the ellipticity and orientation of our 



elliptical fits are i n general consistent with previous! 
publi shed results (IMaughan et al~| | 2008t IDonahue et al 



IDe Filippis et al.l 
IGirardi et alJ 120061: 



iy 
1. 
.1. 



20051 IPiffaretti et al 
. _ McNamara et al.1 1200a 

Schmidt et alJ 120011: iNeumann and Arnaudl 12000 ) . Fi- 
nally, the best-fit concentration parameters from our 
Nagai model fits to three of the five clust ers are consis- 
tent with those found from X-ray data ([Arnaud et al.1 
[2001 ; MS 0015.9+1609 and MS 1054.4-0321 have 
significantly lower values of csoo- We describe each 
cluster in detail below: 

Abell 697: We find that Abell 697 has a significant el- 
lipticity, and it is not well described by either a 
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TABLE 4 

Beta model fit parameters compared to OVRO /BIMA/ Chandra results 



cluster 


yo (io- 4 ) 




r c (arcsec) 


5RA 


<5dec 




Bolocam OV/BI/C/i 




Bolocam OV/Bl/Ch 







Abell 697 


3.63 


± 


0.31 


9 

z - zy -0.24 


0.607 


49 


± 


6 


43 2 +2,1 
^'^-2.0 


-6.2 


± 


4.2 


-14.0 


± 


4.6 


Abell 1835 


2.95 


± 


0.37 


, lg +0.19 


0.670 


49 


± 


6 


32.4+J- 4 


-8.2 


± 


4.4 


1.5 


± 


5.0 


MS 0015.9+1609 


2.80 


± 


0.31 


o ck+0.15 
z - oo -0.15 


0.744 


69 


± 


10 


42.9+^ 


18.2 


± 


5.2 


5.3 


± 


4.7 


MS 0451.6-0305 


3.05 


± 


0.25 


9 79 +0.15 
z - ' -0.13 


0.795 


18 


± 


6 


36.0ll| 


14.2 


± 


4.6 


7.9 


± 


4.2 


MS 1054.4-0321 


1.92 


± 


0.20 




1.083 


70 


± 


12 


70.5l 6 6 ;| 


-1.0 


± 


4.4 


-1.5 


± 


4.6 



Note. — Table of the best-fit parameters and 1<t uncertainties when we fit a spherical beta model to our 
data using the value of j3 found by B06 using OVRO/BIMA and Chandra data (in contrast to our nominal 
beta model fits with j3 = 0.86). From left to right the columns give our best-fit values of yo, the B06 best-fit 
values of yo, the value of /3, our best-fit values of r c , the B06 best-fit values of r c , and the RA and dec offsets 
of the centroids of our best-fit models compared to the X-ray centroids of the B06 fits. Compared to B06, 
we find a significantly larger value of yo for Abell 697, and we find significantly larger values of r c for Abell 
1835 and MS 0015.9+1609. We also find small, but measurable, centroid offsets for all of the clusters other 
than MS 1054.4-0321. 



beta model or a Nagai model. The poor model 
fits result from the cluster appearing significantly 
extended in the SW direction, and extremely com- 
pact in the NE direction. Abell 697 is the only 
cluster in our sample that is not adequately de- 
scribed by the models. The ellipticity we find for 
Abell 697, e = 0.37 ± 0.05, is roughly consistent 
with the ellipticity of — 0-25 found from X-ray 
data (iDe Filippis et al.ll2005l iMaughan et al.ll2008l; 
iGirardi et al.l l2006f K We find a position angle of 
—24 deg for Abell 697, w hich is similar to th e 
value of —16 deg found by IGirardi et al.l ((2006), 
but somewhat misaligned to the p osition angle of 
16 deg found by IDe Filippis et al.l p005h . 

Abell 1835: We detect an ellipticity in our image 
of Abell 1835, and we find that it is well de- 
scribed by either an elliptical beta or Nagai 
model. Our best-fit spherical Nagai model pa- 
rameters with A = 0.9, B = 5.0, and C = 



0.4 (po 



11.1 ± 2.6 x 10~ n erg/cm J 



4.6 ± 1.1 arc min) are consistent with the val- 
ues found in IMroczkowski et aT] (|2009D (po = 
13.6 x 10 -11 erg/cm 3 , r s = 4.3 arcmin) using 
a combination of SZA SZ data and Chandra X- 
ray data. X-ray measurements find an elliptic- 
ity of 0.1 - 0.2 for Abell 1835 (IDe Filippis et~aTl 
[20051 : IMcNamara et al.|[2006l: iSchmidt et al.ll200ii r 
consistent with the value of e = 0.27 ± 0.07 
we find with Bolocam. We find a position an- 
gle of —16 deg for Abell 697, which is similar 
to the values of 7, —2 0, and —30 deg found b y 
IDe Filippis eTaLl ((20051) lMcNamara"etaII (|2006l) . 
and ISchmidt et al.l (|200ll) . 



MS 0015.9+1609: MS 0015.9+1609 appears to be el- 
liptical in our image, and it is well described 
by either a spherical or elliptical model. As 
with Abell 697 and Abell 1835, X-ray data fa- 
vor slightly lower ellipticities, e < 0.20, com- 
pared _to_wbM_wB_find with Bolocam, e = 0.24 ± 
0.08 (IDe Filippis eTaLl 12001 IMaughan et al.1120081: 
Piffa retti et al.l 120031 ). We find a position angle of 
68 deg for MS 0015. 9+1609, fairly close to the value 
of 47 deg found by iPiffaretti et all (|2003l) . but al- 
most orthogonal to the value of —49 deg found by 



IDe Filippis eTaLl (|2005T ). 

MS 0451.6-0305: MS 0451.6-0305 also appears to be 
elliptical in our image and is adequately described 
by either an elliptical beta or Nagai model. We 
find an ellipticity of e = 0.26 ± 0.06 for MS 0451.6- 
0305, in excellent agreeme nt with ellipticities de- 
tcrmincd using X-ray data (|De Filippis et al.ll2005l: 
iDonahue et al.ll2003fi . Additionally, our best-fit po- 
sition angle of 85 deg agrees well w i th the value of 
84 deg found by IDe Filippis et al.l (120051) and t he 
value of —75 deg found bv lDonahue "et&U ((20031) . 

MS 1054.4-0321: Our image of MS 1054.4-0321 shows 
no evidence for ellipticity, and it is well described 
by either an elliptical or spherical model. Al- 
though MS 1054.4-321 does not appear to be el- 
liptical in our data, X-ray data show a clear 
ellipticity oriented along the east-west directio n 
(|Neumann and Arnaudll2000HJeltema et al.ll2001l) . 
However, our non-detection of an ellipticity, e = 
0.09 ±0.07, is only marginal ly inconsistent with the 
X-ray value determined by INeumann and Arnaudl 
(l2000l) . e = 0.29. 

SDS1: We have attempted to fit cluster models to the 
SDS1 map, but we find best-fit amplitudes that are 
consistent with and best-fit scale radii that are 
large compared to the size of our images (i.e., scale 
radii that are large enough to produce profiles that 
are approximately constant over the entire image). 

5. MODEL-INDEPENDENT IMAGES AND Y S z ESTIMATES 

Rather than using models, which at best provide an ad- 
equate description of what clusters look like on average, 
we have chosen to derive observable quantities from our 
images in a quasi-model- independent way. In Section [3. 31 
we described how we calculate a unique transfer function 
for each cluster to quantify the effects of our noise filter- 
ing. In the Fourier space of our images, these transfer 
functions are a set of two-dimensional complex numbers 
that describe how an input cluster image is filtered as a 
function of two-dimensional Fourier mode. We obtain an 
unfiltered, or deconvolved, image of the cluster by Fourier 
transforming our image, dividing by the two-dimensional 
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complex transfer function, and then Fourier transform- 
ing the result back to image spac^H- At the largest scales 
in our map, the transfer function has a magnitude of 
~ 0.2 (see Figure |3]), resulting in a numerically stable, 
but significant, amplification of the large scale noise. In 
particular, the residual atmospheric noise and primary 
CMB fluctuations, which had been filtered to be approx- 
imately white, produce a significant low-frequency noise 
component in our deconvolved images. We estimate the 
noise in our deconvolved images by deconvolving each of 
the 1000 noise realizations for each cluster. Since the 
off-diagonal elements of the noise covariance matrix are 
significant due to the low-frequency noise, we estimate 
all of our measurement uncertainties from the standard 
deviation of measuring the same quantity in each of our 
1000 deconvolved noise realizations (see below and Sec- 
tioning]). Note that most of the measurement uncertain- 
ties in our processed (i.e., filtered) images are computed 
in the same way, even though their noise covariance ma- 
trix is approximately diagonal. 

By deconvolving the two-dimensional transfer function 
of our data processing we obtain unbiased cluster images, 
modulo smoothing with our PSF and the unconstrained 
DC signal level. Although the transfer function was com- 
puted using the best-fit-elliptical Nagai model, and is 
therefore somewhat model dependent, the dependence is 
negligible other than the determination of the DC signal 
offset of the image (see Section 16.11) . We do not attempt 
to recover the information lost due to smoothing by our 
PSF, but we use the value of ST found from our ellipti- 
cal Nagai model fits to restore the correct DC signal level 
to our images. Due to uncertainties in the model itself, 
especially at large radii where there is little or no observa- 
tional data, along with cluster-to-cluster deviations from 
the model, this does introduce a non-negligible model- 
based bias in our images (the typical model uncertainty 
in the DC signal offset is 5 — 10 ijKcmb, see Section |6~T) . 
However, we emphasize that this bias only affects the DC 
signal level of the images; the shapes of the cluster pro- 
files are essentially model-independent. 

We have computed a model-independent value for Y500 
from these deconvolved images, which is the integrated 
y within rsoo (i.e., the cylindrical F500 rather than the 
spherical Y500, see Table [5]). As mentioned above, since 
our assumption that the noise covariance matrix is di- 
agonal fails for the deconvolved images, we estimate the 
uncertainty in our estimate of Ysz via the scatter among 
the Ysz values determined from each of our noise real- 
izations. As an example of the amount of low-frequency 
noise present in our deconvolved images on the scale of 
rsoo, note that the measurement uncertainty on Y500 is 
approximately 10 times larger than it would be if the 
noise was white. For the five clusters in our sample, we 
are able to determine Ysz with an uncertainty of ~ 10%, 
limited mainly by systematics in determining the DC sig- 

7 This method can be compared to the deconvolution method 
emplo yed by APEX-SZ for their analysis of Abell 2163 and Abell 
2204 INord et al.ll2009l ; iBasu et al.ll2010ft . They first determine 
what a point-like object looks like in their image after being filtered. 
Next, they fit this filtered point-source image to the map pixel 
with the largest S/N and subtract it from the image. The process 
is repeated until the map is consistent with noise; the sum of all 
the unfiltered point-like images removed from the map gives the 
deconvolved cluster image. 



nal offset and flux calibration. 

We examine the self-similar scaling that is expected 
between Ysz and the total cluste r mass using the relation 
given m Bon amente et all (|2008[) . 

Y S zD\E{zr^ cx f-*J 3 M 5 g ll 

where Da is the angular diameter distance, E{z) = 

yj£l m (l + z) 3 + Qa, and we have as sumed M nns is a 
good proxy for the total cluster mass (|Allen et al.lkuuv 
iMantz et al.ll2010al ). T he scatter in this rel ation is ex- 
pected to be < 10% (|Kravtsov et al.ll2006h . and cur- 
rent meas urements are rough l y consistent with this 
prediction dMorandi et aT1l2007t iBonamente et al] 120081: 
Marrone et alj l2009t" Plagge et all 120101: iHuang et at] 
2010t lAndersson et all 12010ft . Given our small sample, 
we do not attempt to constrain the intrinsic scatter in 
the Y — M r e lation , but we follow the formalism of 
iMarrone et al.l (|2009f) and P10 to fit a logarithmic scal- 
ing of the form y = a + bX with the intrinsic scatter set 
to 10%. For the Y500 — M gas 500 relation we find best-fit 
values of a = -5.46 ± 0.84 and b = 1.63 ± 0.71, con- 
sistent with the self-similar prediction of b = 5/3. The 
overall scatter of our data about the fit is 13%, consis- 
tent with an intrinsic scatter of ~ 10% given our ~ 10% 
uncertainty on Ysz (see Figure H]). 

Additionally, our results are consistent with Ysz ~ 
M gas scaling relations measured by other groups. For 
example, our results using model-independent Ysz esti- 
mates from Bolocam and M gas estimates from Chandra 
within rsoo agree well with the results in P10 using Ysz 
estimates from the SPT best-fit beta model and M gas es- 
timates primarily from XMM-Newton (a = —5.73 ±0.43 
and b = 2.12 ± 0.45 within r 50Q and a = -5.92 ± 0.41 
and b = 1.97 ± 0.44 within r 25 oo)- The results in 
IBonamente et all (|2008[ L using beta model fits to 30 GHz 
OVRO/BIMA SZ data and Chandra X-ray data within 
r25oo, also match our results quite well (a = — 5.22 ± 1.77 
and 6= 1.41 ±0.13). 

6. TESTS FOR SYSTEMATIC ERRORS 

We are primarily concerned with two types of system- 
atic errors that may occur in our analysis: those caused 
by our (minimal) use of a cluster model, specifically an 
elliptical Nagai model, and those caused by characteris- 
tics of our noise that are not properly accounted for in 
our analysis. We have run extensive tests to quantify the 
level of systematic error we can expect from each of these 
two sources. The details of these tests are described be- 
low. In the end, we find that the amount of systematic 
error in our images is negligible, with the exception of the 
estimation of the DC signal offset in our images using the 
Nagai model. 

6.1. Model dependence of results 

Since we compute the transfer function for each cluster 
using the best-fit elliptical Nagai model for that cluster, 
our deconvolved images necessarily have some model de- 
pendence. In order to quantify the amount of model 
dependence, we have computed transfer functions for a 
range of elliptical Nagai models for one of the clusters 
in our sample, MS 0451.6-0305. Relative to the best-fit 
model, we have varied the scale radius, r s , the ellipticity, 
e, the position angle, and the centroid location, SKA 
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TABLE 5 
Model-Independent Y500 estimates 







JK noise 


CMB/PS 


map DC signal 


flux cal. 


cluster 


*500 


cry 


cry 


cry 


cry 


Abell 697 


59.2 ±6.3 


1.2 


0.5 


5.6 


2.5 


Abell 1835 


47.5 ±4.9 


1.3 


0.5 


4.3 


2.0 


MS 0015.9+1609 


23.3 ±2.2 


1.6 


0.6 


0.9 


1.0 


MS 0451.6-0305 


24.7 ±2.0 


1.2 


0.6 


0.9 


1.1 


MS 1054.4-0321 


8.8 ± 1.0 


0.8 


0.4 


0.2 


0.4 



Note. — Model independent estimates of Y50O1 in units of 10 — 1 ster, along 
with the associated uncertainties. From left to right the columns give the value and 
total uncertainty of our estimate of I50O1 the uncertainty due to our jackknife noise 
model (i.e., excluding CMB and point sources), the uncertainty due to CM B and 
point sources, the uncertainty in the DC signal level of the map (see Section 16.111 . 
and the uncertainty in our flux calibration. Note that the jackknife noise estimate 
is dominated by large angular-scale residual atmospheric fluctuations; the noise on 
beam-size scales is negligible. 
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Fig. 6. — Scaling relation between Y$z and M gas for the five 
clusters in our sample. The black squares give the Bolocam model- 
independent cylindrical Y500 estimates and Chandra M ga s,500 es- 
timates (see Table [TJ, and the black dashed line represents the 
best fit to these data. We also show fits of a similar type obtained 
by various authors: red dashed: SPT beta-model-derived Y500 vs. 
Mgas, 500i primarily obtained from XMM-Newton data (P10); solid 
blue: same, but within r25oo; dot-dashed green: OVRO/BIMA 
beta- model-derived V2500 vs. M ga s,2500 obtained from Chandra 
data (Bonamcntc ct al. 2008). The difference between the red- 
dashed line and our data is likely caused by systematic differences 
between model-derived and model-independent estimates of Ysocfl. 
The scatter of our data relative to our best fit is 13%, consistent 
with the expected intrinsic scatter of ~ 10% given our ~ 10% un- 
certainty on Y5oo- 

a P10 calculated both model-derived and model-independent es- 
timates of Y§z f° r the 15 clusters in their sample. Within a radius 
of ^2500 1 the two estimates of Y2500 are on average consistent with 
each other. However, the P10 model-derived estimates of Y500 are 
on average a factor of 1.5 larger than the model-independent esti- 
mates for the 10 clusters used to determine the Ygz — Mgas scaling 
relation. As a result, the P10 model-derived scaling relation for 
Vgoo — Mgas,500 will be systematically higher by log 10 (1.5) ~ 0.2 
compared to a model-independent scaling relation. 

and (5dec, by increasing and decreasing each parameter 
individually by its la uncertainty. We then deconvolved 
our processed map of MS 0451.6-0305 using the transfer 
function computed from each model and subtracted the 
resulting map from the one produced using the transfer 

8 The power law slopes (A., B, and C) were held fixed for all of 
our model fits. Due to the large degeneracy between these values 
and r s , we have effectively included variations in the power law 
slopes by varying the value of r a . 



function for the best-fit model. In each case, the residual 
map was approximately white, with an RMS of 1.5, 0.6, 
and 0.5 ^Kcmb for variations in r s and e, 9, and SKA 
and Sdec, respectively. Since the typical noise RMS of 
our deconvolved maps is ~ 10 [iKcmb, the additional 
RMS introduced by our uncertainty in determining the 
model used for calculating a transfer function is quasi- 
negligible. Note that the best-fit elliptical Nagai model 
will not provide an exact description of a real cluster. 
However, the elliptical Nagai model does provide an ad- 
equate description of four of the five clusters we have 
observed, indicating that the difference between the true 
cluster profile and the model profile is in general less 
than our noise. Therefore, the artifacts in our decon- 
volved map produced by using a model to describe the 
cluster will be smaller than the artifacts produced by our 
measurement uncertainty on the best-fit model. 

Additionally, we created a deconvolved map of MS 
0451.6-0305 using the transfer function for a point-like 
source. The resulting profile is significantly different 
from the profile obtained using the transfer function for 
the best-fit Nagai model, indicating that the naive cal- 
culation of a transfer function using a point source is 
inadequate. Compared to using the transfer function 
calculated from the best-fit elliptical Nagai model, the 
peak decrement is reduced by ~ 50 /jKcmb, while the 
magnitude of the SZ signal at the edge of the map is 
increased by ~ 50 (jKcmb (i.e., the deconvolved cluster 
image from a point-source transfer function is systemat- 
ically broader). 

As described in Section |H we use the best-fit el- 
liptical Nagai model to determine the DC signal level 
in our maps since its value is unconstrained by our 
data. The measurement uncertainty in the value of 
the DC signal is small, < 1 ^Kcmb, but there is 
a large amount of uncertainty i n the model at larg e 
radii. Ba sed on the results from iBorgani et al.l (|2004l) . 
N07, and IPiffaretti and Valdarninil (|2008D presented in 
lArnaud et all 1)20091 ) . the RMS scatter in the pressure 
profiles from cluster to cluster in simulations is < 25%. 
Therefore, we include an additional 25% systematic un- 
certainty on the DC signal that we add to the decon- 
volved map. Specifically, we estimate that the model 
uncertainty in the DC signal level of our map is 25% of 
the signal level at the edge of the map. 
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6.2. Noise characteristics 

As mentioned in Section [31 we make the approximation 
that the noise covariance matrix is diagonal in our pro- 
cessed maps (i.e., there are no noise correlations between 
pixels). Although fluctuations in both the atmospheric 
emission and the CMB are correlated over many pixels, 
the high-pass filter we apply to our data timestreams 
eliminates these correlations within our ability to mea- 
sure them. To test for noise correlations, we added pro- 
cessed cluster images of the best-fit models to the 1000 
noise realizations for each cluster (i.e., for a given clus- 
ter we separately added each of the four best-fit cluster 
models from Tables [5] and [3] to each of the 1000 noise re- 
alizations). We then fit a model of the same type (e.g., if 
we added the best-fit elliptical Nagai model to the noise 
realization, then we fit an elliptical Nagai model to the 
cluster-model-plus-noise map) to each of these model- 
cluster-plus-noise maps and examined the distribution of 
X 2 values for these fits. For all five clusters, we found that 
the distribution of x 2 values matched the predicted x 2 
distribution obtained from the diagonal noise covariance 
matris0. As an example, the fit quality of our measured 
X 2 distribution for the elliptical-Nagai-model-plus-noise 
realizations to the predicted x 2 distribution, quantified 
by a PTE, is 0.58, 0.02, 0.08, 0.66, and 0.21 for Abell 
697, Abell 1835, MS 0015.9+1609, MS 0451.6-0305, and 
MS 1054.4+321. (see Figured]). Therefore, we conclude 
that there is a negligible amount of correlated noise in 
our processed images and our assumption that the noise 
covariance matrix is diagonal is valid. Furthermore, we 
have estimated all of our parameter uncertainties (po, r s , 
etc.) directly from the distribution of values calculated 
from our noise realizations. Consequently, any failure of 
our assumption that the covariance matrix is diagonal for 
the processed images will only affect the pixel-weighting 
and x 2 values for our model fits. 

However, when we perform the same test of fitting a 
model to our model-cluster-plus-noise maps using noise 
realizations that have been deconvolved with our trans- 
fer function, the result is significantly different. The dis- 
tribution of x 2 values calculated from our deconvolved 
noise realizations is significantly broader than the pre- 
dicted distribution based on uncorrelated noise (see Fig- 
ure [?])• This result is not surprising since the decon- 
volution enhances the large-scale signals in the images, 
including residual atmospheric noise and CMB fluctua- 
tions. Since the noise in the deconvolved images is sig- 
nificantly spatially correlated, our assumption that the 
noise covariance matrix is diagonal fails. Therefore, we 
estimate the uncertainties for the deconvolved images us- 
ing the spread in values for the noise realizations rather 
than from the diagonal elements of the noise covariance 
matrix (e.g., the uncertainties in the radial profiles are 
determined from the RMS spread in the radial profiles of 
th e noise realizations ). T his is the same tec hnique used 
by iNord et al.1 (|2009f ) and iBasu et al.1 (|2010t l to analyze 
APEX-SZ data. 

The model fits to the model-cluster-plus-noise realiza- 

9 We have fit for N params free parameters in our fits to both 
the actual data and the noise realizations, with N para vns — 5 or 
7 for the spherical or elliptical fits. Therefore, the predicted \ 2 
distribution is for N p i x — N pararns DOF, where N p i x is the number 
of map pixels. 



tions also provide us with estimates of the uncertainties 
and biases associated with our model-parameter fitting. 
Specifically, we obtain 1000 best-fit values for each pa- 
rameter; the standard deviation of these values then gives 
the uncertainty on our estimate of that parameter in our 
actual data map. These uncertainties are given in Ta- 
bles [5] and [3J In addition, if our parameter estimation 
algorithm is free from biases, then we should, on aver- 
age, recover the parameters of the input model that was 
added to the noise realizations. In practice, we find a 
small, but measurable, bias in our estimates of the pres- 
sure normalization and scale radius in our fits; the bias 
is typically < 10% of the uncertainty on each paramter. 
We find no measurable bias in our estimates of the other 
fit parameters. 

Additionally, we have used our signal- free SDS1 maps 
to further verify our model-fitting proceedure and to 
search for any components of the noise that have not 
been included in our noise estimate. First, we inserted 
model clusters into the SDS1 data timestreams based on 
the best-fit elliptical Nagai profile for each of the five clus- 
ters in our sample. These data were then processed and 
an elliptical Nagai model was fit to each resulting image. 
In each fit, there are 6 free parameters (po, r s , e, 9, <5RA, 
and <5dec), giving us a total of 30 fit parameters for the 
5 model clusters. Of these 30 fit parameters, 17 (57%) 
are within la of the input value, 26 (87%) are within 2a 
of the input value, and all 30 are within 3cr of the input 
value; these results indicate that our model fitting and 
parameter error estimation are working properly. Addi- 
tionally, we obtain a reasonable goodness of fit for the 
models using these best-fit parameters, quantified by a 
PTE ~ 0.8, providing further evidence that there is no 
significant noise in the data that has not been included 
in our noise estimate. Note that since the five cluster 
model profiles are fairly similar, we obtain comparable 
PTEs for all five profiles. 

7. SUMMARY 

We have presented the first results from our program 
to image the SZ effect in galaxy clusters with Bolo- 
cam. These images have a beam-smoothed RMS of 
~ 10 /jKcmb, and a resolution of 58 arcsec. Given this 
noise level, we are able to measure SZ signal in radial 
profiles to approximately the edge of our maps, which 
corresponds to 6-7 arcmin or 1 — 2 times ^oo- In order 
to subtract noise from atmospheric fluctuations, we ef- 
fectively high-pass filter our cluster images. However, we 
are able to deconvolve the effects of this filter with bi- 
ases that are negligible compared to our noise level, other 
than our recovery of the DC signal level. In fitting our 
images to spherical and elliptical beta and Nagai models, 
we find no preference between the beta and Nagai mod- 
els due to the degeneracy between these models over the 
angular range to which our data are sensitive, but our 
data do show a definitive preference for elliptical models 
over spherical models. The weighted mean ellipticity of 
the five clusters is e = 0.27 + 0.03, consistent with results 
from X-ray data. Additionally, the best-fit model pa- 
rameters we determine from our data are consistent with 
those found from previous X-ray and SZ measurements. 
We have also obtained model-independent estimates of 
Ygzi arL d we find scaling relations between Ysz and clus- 
ter mass that are consistent with self-similar predictions, 
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Fig. 7. — Histograms of the \ 2 value for 1000 separate noise realizations for MS 0451.6-0305, overplotted in green with the predicted 
distribution assuming the noise covariance matrix is diagonal (i.e., there is no correlated noise). For each noise realization the best-fit 
elliptical Nagai model profile for MS 0451.6-0305 is added to the noise realization, an elliptical Nagai model is fit to this model-cluster- plus- 
noise realization, and the value of \ 2 ls computed based on the assumption that the noise covariance matrix is diagonal. The vertical red 
line shows the value of \ 2 f° r the actual data for MS 0451.6-0305. The left histogram shows the processed data, and the right histogram 
shows the data for the deconvolved image. The predicted and actual \ 2 distributions for the processed data overlap, indicating that there 
are minimal correlations between map pixels. However, the actual x 2 distribution for the deconvolved image data is much broader than 
the predicted distribution, indicating that there are significant noise correlations between map pixels in the deconvolved images. 



with a scatter that is consistent with expectations for a 
~ 10% intrinsic scatter. 
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APPENDIX 



This appendix includes images and radial profiles of the processed and deconvolved maps, an image of the processed 
residual map after subtracting the best-fit elliptical Nagai model, and an image of one of the 1000 noise estimates 
generated via jackknife realizations of the data and a model for the astronomical noise. We have smoothed all of the 
images using a Gaussian beam with a FWHM of 58 arcsec. A white dot representing the FWHM of the effective PSF 
for these beam-smoothed images is given in the lower left of each image. The solid white contour lines in the images 
represent a S/N of —2, —4, .., and the dashed white contour lines represent a S/N of +2, +4, .... The deconvolved images 
contain a significant amount of noise that is correlated over large angular scales, along with a model-dependent DC 
signal offset, and we therefore do not display noise contours on the deconvolved images. The error bars on the radial 
profiles are estimated from the spread in radial profiles computed from our noise realizations, and therefore do include 
all of the large-angular-scale noise correlations (although they do not include the uncertainty in the DC signal level of 
the image). Note that the radial profile bins for the deconvolved images are correlated due to the large- angular-scale 
noise present in those images. 
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Fig. Al. — Abell 697; from left to right and top to bottom we show the deconvovled image of the cluster, the processed image of the 
cluster, the residual map between the processed image of the cluster and the best-fit elliptical Nagai model, one of the 1000 noise realizations 
for the processed data, and a binned radial profile. The contour lines represent a S/N of 2, 4, ... 
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Fig. A2. — Abell 1835; from left to right and top to bottom we show the deconvovled image of the cluster, the processed image of the 
cluster, the residual map between the processed image of the cluster and the best-fit elliptical Nagai model, one of the 1000 noise realizations 
for the processed data, and a binned radial profile. The contour lines represent a S/N of 2, 4, ... 
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Fig. A3. — MS 0015.9+1609; from left to right and top to bottom we show the deconvovled image of the cluster, the processed image 
of the cluster, the residual map between the processed image of the cluster and the best-fit elliptical Nagai model, one of the 1000 noise 
realizations for the processed data, and a binned radial profile. The contour lines represent a S/N of 2, 4, ... 
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Fig. A4. — MS 0451.6-0305; from left to right and top to bottom we show the deconvovled image of the cluster, the processed image 
of the cluster, the residual map between the processed image of the cluster and the best-fit elliptical Nagai model, one of the 1000 noise 
realizations for the processed data, and a binned radial profile. The contour lines represent a S/N of 2, 4, ... 
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Fig. A5. — MS 1054.4-0321; from left to right and top to bottom we show the deconvovled image of the cluster, the processed image 
of the cluster, the residual map between the processed image of the cluster and the best-fit elliptical Nagai model, one of the 1000 noise 
realizations for the processed data, and a binned radial profile. The contour lines represent a S/N of 2, 4, ... 
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Fig. A6. — SDS1; from left to right and top to bottom we show the processed image of the field, one of the 1000 noise realizations for 
the processed data, and a binned radial profile. The contour lines represent a S/N of 2,4,... The thin grey lines show the radial profiles 
for each of the 1000 noise realizations. 



